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ABSTRACT 

Newly-born pulsars offer favorable sites for the injection of heavy nuclei, and for their further 
acceleration to ultrahigh energies. Once accelerated in the pulsar wind, nuclei have to escape from the 
surrounding supernova envelope. We examine this escape analytically and numerically, and discuss 
the pulsar source scenario in light of the latest ultrahigh energy cosmic ray (UHECR) data. Our 
calculations show that, at early times, when protons can be accelerated to energies E > 10^° eV, the 
young supernova shell tends to prevent their escape. In contrast, because of their higher charge, iron- 
peaked nuclei are still accelerated to the highest observed energies at later times, when the envelope has 
become thin enough to allow their escape. Ultrahigh energy iron nuclei escape newly-born pulsars with 
millisecond periods and dipolc magnetic fields of ~ 10^^~ G, embedded in core-collapse supernovae. 
Due to the production of secondary nucleons, the envelope crossing leads to a transition of composition 
from light to heavy elements at a few EcV, as observed by the Auger Observatory. The escape also 
results in a softer spectral slope than that initially injected via unipolar induction, which allows for 
a good fit to the observed UHECR spectrum. We conclude that the acceleration of iron-peaked 
elements in a reasonably small fraction (< 0.01%) of extragalactic rotation-powered young pulsars 
would reproduce satisfactorily the current UHECR data. Possible signatures of this scenario are also 
discussed. 



1. INTRODUCTION 

The origin of the highest energy cosmic rays 
still remains a mystery (see iKotera fc Olintd 120111 : 
iLetessier-Selvon fc Stanevll2011l for recent reviews) . The 
measurement of a fl ux suppression at the hig hest energies 
(jAbbasi et al.ll2008t [Abraham ct al."2010bV reminiscent 
of the "GZK cut-ofP' ([Greiscn 1966; Zatscpin & Kuzmin 
Il966f ) produced by the interaction of particles with the 
cosmic microwave background (CMB) photons for prop- 
agations over intergalactic scales, has appeased the de- 
bate concerning the extragalactic provenance of UHE- 
CRs. This feature not only suggests that UHECRs would 
originate outside our Galaxy, but also that the sources 
of the highest energy particles should be located within 
^ 100 Mpc distance, in our local Universe. However, the 
sources remain a mystery and results from the Auger Ob- 
servatory on the arrival directions and chemical compo- 
sition of UHECRs make the picture even more puzzling. 

Hints of anisotropies in the sky distribution of cosmic 
rays above 60 EeV were reported by the Auger Obser- 
vatory, but most of the anisotropy signal seems to issue 
from a clustering of events over a few tens of degree s 
around the region of Centaurus A (jAbreu ei~aII[20T?l . 
No powerful sources are observed in the direction of the 
highest energy events. This might be explained by strong 
deflections that cosmic rays could experience in presence 
of particularly intense extragalactic magnetic fields or if 
they were heavy nuclei. This absence might also find a 
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natural explanation if the sources were transient, such 
as gamma-ray bursts or newly-born pulsars. The de- 
fiection in the extragalactic magnetic fields should in- 
deed induce important time delays (~ 10^ yr for one 
degree deflection over 100 Mpc) between charged parti- 
cles and the photons propagating in geodesies, so that 
the sources should already be extinguished when cos- 
mic rays are detected on Earth. Even in this case, for 
proton dominated compositions and intergalactic mag- 
netic fields of reasonable strengths, the UHECR arrival 
directions are expected to trace the large scale struc- 
tures where the transient sources are distributed, with a 
possible bias (jKalli ct al. 201l]). The precise role of ex- 
tragalactic magnetic fields in UHECR propagation may 
be clarified in the future through extensiv e Faraday ro- 
tation surveys (see, e.g.. iBeck et al.l[2007f ) and indirect 
measurements of gamma- ray halos around blazars (e.g., 
iNeronov fc Semikod [20091 ). 

The composition measurements at the highest energies 
of the A uger Observatory are surprising. [Abraham et al.l 
(|2QlQa| ) report a trend from a proton dominated compo- 
sition at a few EeV toward an iron dominated compo- 
sition at around 40 EeV (continuing up to 60 EeV, see 
lAbreu et al.ll2011bD . assuming that hadronic interaction 
models can be extrapolated to these energi es. This trend 
is not confirmed by the HiRes experiment (jAbbasi et al.l 
|2^ nor by the preliminary data of the Telescope Array 
(jTanicda ct al. 20111), who report light primaries in the 
Northern hemisphere (while Auger observes the Southern 
hemisphere). One may note however that both results re- 
main consistent with those of Auger within quoted sta- 
tistical and systematic errors. 

From a propagation point of view, heavier nuclei 
are favored compared to light elements for a given en- 
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ergy as they can travel hundreds of megaparsecs be- 
fore losing their energy by photo-disintegration processes 
on the cosmic backgrounds due to their lower energy 
per barvon (e.g..lStecker fc Salamonlll999l:lBertone et alJ 
[20nllAllard et alJ 120051 . 1200a iHooper et al.ll2n05h . Nu- 
clei of charge Z can also be in principle accelerated to 
an energy typically Z times larger than protons in a 
given electromagnetic configuration. Propagation mod- 
els where a heavy composition arises at the highest 
energies due to a combination of a low proton maxi- 
mum acceleration energy (around 10 EeV) and Z times 
higher maximum energies for heavier elements (present 
in a slightly higher abundance than Galactic) have been 
shown to reproduce the composition trends o bserved by 
Auger (jAUard et al.ll2008HAloisio et al.ll2009D . However, 
these works focus on the propagation, and do not pro- 
vide a plausible source for the injection of these specific 
compositions. The problem of finding powerful sources 
that inject mainly these low abundance elements and of 
their escape from the acceleration site remains open. 

Heavy nuclei dominated injection models are quite 
rare in the astrophysical literature of candidate sources. 
A direct injection of large proportions of heavy nu- 
clei into an acceleration region requires either an ini- 
tial metal-rich region, or an efficient nucleosynthesis in 
the accelerating outfiow. These requirem ents are hardly 
met by fireball-ty pe gamma-ray bursts (jLemoind 120021 : 
iPruet et al.]l2002t ). Active galactic nuclei (AGN), which 
are the other popular sites for UHECR acceleration mod- 
els, are observed to have a solar to super-solar metallici- 
ties, but wit h a low proportion of nuclei heavier than n i- 
trogen (e.g., IGroves et al.l [20061: IMathur fc Fields! [2001) . 
Young neutron stars on the other hand possess iron-rich 
surface and early conditions that are propitious for heavy 
nuclei injection. 

Pulsars have been suggested as po ssible accelerators 
of co smic-rays since their discovery (jGunn fc Ostriked 
|1969( ). due to their important rotational and magnetic 
energy reservoirs. Galactic pulsars have been suggested 
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200J . iBlasi et alJ(|2000l ) proposed that iron nuclei accel 



erated in the fastest spinning young neutron stars could 
explain the observed cosmic rays above the ankle in a 
Galactic source scenario. They assumed that the strip- 
ping of heavy nuclei from the surface of the star is a plau- 
sible seeding and derived a spectru m based on th e spin 
down of young pulsars (J oc E~^). lAronsI (|2003| ) stud- 
ied the birth of extragalactic magnetars as the source 
of ultrahigh energy protons, developing the acceleration 
mechanism in detail and assuming that the magnetar 
wind disrupts the supernova envelope to allow the es- 
cape of accelerated particles. 

The iBlasi et al.l ([2000l ) and lArond ([200l proposals 
for the origin of UHECRs were elaborated to explain 
the absence of the GZK cut -off in the observed spec- 
trum reported by AGASA (jTakeda et al.l I1998D with- 
out invoking the so-called top-down models (see, e.g., 
iBhattachariee fc Sigll [2000l ) . An increase in the expo- 
sure at the ultrahigh energies by the HiRes and Auger 
Observatories have shown that the UHECR spectrum 



is consistent with a GZK cutoff (jAbbasi et al.l 120081 : 
[Abraham et al.|[2"010bD . A decade ago, the chemical com- 
position was also barely detectable at the highest ener- 
gies while recent results suggest a puzzling trend toward 
heavier nuclei. A new investigation of the pulsar scenario 
as UHECR sources is timely, in the light of the data that 
has been recently acquired. 

In this paper, we examine the key mechanisms involved 
in the production of UHECRs by newly-born pulsars, 
and discuss their implications, considering the latest ob- 
servational results. We focus in particular on the effects 
of the dense supernova envelope that surrounds the neu- 
tron star, and that accelerated particles have to traverse 
on their way to the interstellar medium. We perform 
detailed analytical and numerical Monte-Carlo calcula- 
tions of the envelope crossing and predict the out-coming 
features that particles should bear after the escape. It 
is found that a small fraction of extragalactic rotation- 
powered young pulsars embedded in supernova: could 
satisfactorily explain the latest UHECR observations. 

The layout of this paper is the following. In Section [21 
we review and update the discussions related to the pro- 
duction of UHE heavy nuclei in newly-born pulsars. In 
Section [3l we describe the supernova envelope modeling 
used to develop our analytical estimates and to perform 
our numerical simulations of the escape of UHECRs. 
Our main results on the escape of UHECRs from the 
supernovffi envelopes are presented in Section [3] In Sec- 
tion 2] we discuss the implications of the newly-born pul- 
sar model in view of the available UHECR observations. 
There, we argue how a reasonably small fraction of extra- 
galactic fast spinning young pulsars embedded in super- 
novae could reproduce satisfactorily the current UHECR 
data, and discuss observable signatures that could probe 
the pulsar model. Our conclusions are drawn in Sec- 
tion [5] 

2. UHE HEAVY NUCLEI PRODUCTION 
IN NEWLY-BORN PULSARS 

In this section, we review and discuss some key points 
related to the production of UHE heavy nuclei in newly- 
born fast-spinning neutron stars. Our numerical applica- 
tions focus on isolated rotation-powered pulsars of radius 
i?*_io = i?,/10km, angular velocity Vl^ = il/W^s^^, 
principal moment of inertia J45 = I/IO^^ gcm^, and 
magnetic dipole moment ^^30.5 = /i/lO'^"'^ cgs with /i = 
BRl/2 = 1030-5 cgs(B/6 X lO^^ G)Rlj^„, with B the 
surface dipole field strength. We show in Section |3] that 
such parameters would enable the escape of UHE nuclei 
from the surrounding supernova envelope. 



2.1. Acceleration by unipolar induction 

Rapidly rotating neutron star magnetospheres 
are promising particle a, cceleration sites (see, e.g., 
iShapiro fc Teukolskvl 119831 and references therein). In 
the out-flowing relativistic plasma, the combination 
of the fast star rotation and its strong magnetic field 
can induce, in principle, potential differences of order 
$ = rt^fi/c^. Provided that particles of charge Z can 
experience a fraction 77 of that potential, they can 
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be ac celerated to the energy (jBlasi et alj 120001 : lAronsI 
[2OOI : 

E{n) = Ze $ 77 = 3 X lO^o Z26 m M30.5 cV (1) 
where rji = rj/OA and Z26 = Z/26 for hon nuclei. 

Energy losses by gravitational waves and electromag- 
netic radiation lead to the spin-down of the pulsar (see 
iShapiro fc Teukolskvlll983l and references therein j3, and 
thus to the production of particles of lower and lower 
energies as time goes. Under th e assumption that 
the G oldreich- Julian charge density (jGoldreich fc JulianI 
Il969f l is entirely tapped in the outflow for acceleration, 
and using the expression of the pulsar spin-down rate, 
one can derive the en ergy spectrum of the accelerated 
particles (|Aronsll2003l) : 
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where Eg is the critical gravitational energy at which 
gravitational wave and electromagnetic losses are equal. 
The gravitational wave losses start dominating at the 
highest energies when the magnetic field of the star be- 
comes ^ > 10^'^ cgs. Magnetars are thus affected by these 
losses. For pulsars with milder fields that are the main 
concern of this paper, gravitational wave losses are negli- 
gible, and - Eg ^ 10^° eV . In this case, the injected spec- 
trum reads (jBlasi et al.ll2000t) : 

^ = 5 X 10^3 /45(Z26 M30.5 E^^y^cY-', (3) 
dE 

The spin-down time at which particles of energy E can 
be accelerated in the voltage dr op, when gr avitational 
wave losses are neglig ible, reads ()Aron£ll2003[ ): 
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where Ei is the maximum acceleration energy corre- 
sponding to the initial angular velocity Hi. The spin- 
down time at which particles of energy E can be accel- 
erated does not depend on the initial rotation velocity of 
the neutron star fli, for E <^ E-^. 

2.2. Acceleration sites 

Various authors have discussed particle acceleration in- 
side the light cy linder of pulsars and magnetars (see, e.g., 
iHarding fc Lail[2006l for a review). Possible sites include 
the polar ca p region, just above the magnetic pole of 
the s tar (e.g.. lSturro"cMll971l : lHarding fc Muslimovll200lL 
|2002| ). the "slot gap" region along the last ope n field 
line b etween the polar cap and the light cylinder ([Aronsl 
|1983() . an d in the outer gap reg io n close to the light cylin- 
der (e.g.. iCheng et al.l Il986al lbl: iBednarek &: Protherod 
[19971 12002h . Energy losses by curvature radiation are 

^ Numerical simulations of magnetized neutron star relativistic 
winds suggest that the spin-down rate m ay be faster than obtaine d 
in the standard "vacuum dipole" model HBucciantini et aI1l2006l) . 



however likely to prevent the acceleration of particles 
to the high est energies both i n the pola r cap and the 
outer gap. iVenkatesan et al.l (|1997[ ) and lAronsI (|2003[ ) 
discussed that particles accelerated in the wind region 
with r ^ Ri^ with _Rl the radius of the light cylinder, do 
not suffer curvature radiative losses. 

In t h e nex t paragraphs, we follow the arguments of 
lAronsI ()2003[ ) to calculate the radius at which particle 
acceleration is most likely to occur. We also take into 
account the effects of curvature radiation of pions that 
was not previously considered, though it could be more 
constraining than the curvature radiation of photons. 

Outside the light cylinder, the dipole field struc- 
ture cannot be causally maintained and the field be- 
comes mostly azimu thal, with field lines spiraling out- 
wards (jMichell [1991 . In regions of the wind where 
the rest mass density is not dominated by electron and 
positron pairs, the plasma can be considered as force- 
free. In such regions, and for th e case of aligned rotators, 
IContopoulos fc KazanasI (j2002D calculated that charged 
particles flow out with a motion along the (nearly az- 
imuthal) magnetic field lines that becomes negligible 
when r 3> J'min.iin = 7l-Rl- The intial Lorentz factor 
of the particles entering the wind, 7l, can take values 
between 10 — lO'^ depending on the magnetospheric pa- 
rameters. Beyond r ^ ^'min.iin, particles flow out nearly 
radially (they "surf-ride" the fields) and the wind acts 
like a linear accelerator: the Lorentz factor of the out- 
flowing plasma increases linearly as 7w ^ r/Ri^. 

lAronsI (I2003D e xt ended the work of 
IContopoulos fc KazanasI (|2002| ) to oblique rotators 
and to regions in the wind where magnetic dissi- 
pation occurs (i.e., in non force- free regimes), for 
'T > ''diss ^ 2K±i?L- Here k± is the ratio between the 
number density of heavy ions (that we assume equal to 
the Goldreich- Julian density) and of electron-positron 
pairs. Calculations of pair creation i n ultra-magnetized 
neutr oi i stars s ugges t k± ^ 10 — 100 ([Baring fc Harding 
l2001f ). I Aronsl ([20031 ) discussed that surf-riding acceler- 
ation can still occur in these more general cases. He 
argues further that magnetic dissipation via Alfven 
wave emission beyond rjiss would lead to an even more 
efficient surf-riding process, the waves acting as strong 
pondcrmotive forces on the ions. The Lorentz factor of 
the ions (of mass mj) would then reach values as high 
as 7i = Ze77$/(mic2) > 7^ for r > rdiss- The results 
obtained for the unipolar induction toy-model described 
in Section [2.11 can then be applied. 

The curvature rad ius of a surf -riding ion at distance 
r > rmin.iin reads ([Arond [2003[): Pc = 2^17^, where 
p\ ^ rjr is the Larmor radius of the particle[3 One 
can calculate that, to avoid photon curvature radiation 



^ T h e com plete expression of the curvature radiation given by 
I Aronsl II2003I ) is pc = 2pi7^/ cos(f2, /^). The angle between the 
rotation axis and the magnetic dipole moment needs to satisfy 
(f2,/^) < 90° to avoid curvature radiations. In such a configu- 
ration, one can expect an outflow of ions to form from the po- 
lar cap to the rotational equato r, along the last clo sed fie ld lines 
(the so-called "return current " , IGoldreich fc Juli an 1969; Michej) 
[19751 : [Contopoulos etldllTggSl V In the model of.Arons (200:3), it 
is specifically this current of ions that is tapped into the wind for 
acceleration. 
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losses, the acceleration of particles at E21 = i?/10^^eV 
needs to take place at radius greater than: 
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^6.m'E',CzU'A-,'^',;'/'n-'/\rn. (7) 

The cooling timescale for cu rvature radiation of pio ns 
is more constraining; it reads (jHerpav fc Patk6"sll2008f ): 
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where x = E'^h/ {pcA^m^c'). We present here only the 
case of charged pions 7r+ , as this p rocess dominates th e 
case of the emission of TT~ and TfO ()Herpav et alj [20081 ) . 

One can readily see that X - 13£;|iA^6^77j-^^]4(i?LA)^ < 
1 and thus, ^c.tt 3> iacc, for sufficiently large r ^ Ri^ in 
the wind. Numerically, for the same parameters as in 
Eq. ([H]), the acceleration above £'21 = i^/lO^^eV needs 
to take place at r > rmin.c.-n- 2 x 10^ cm to avoid energy 
losses through curvature radiation of charged pions. 

The radiation fields in the pulsar wind are unlikely to 
impact the acceleration of UHECRs. The early neutrino- 
driven wind should end wit hin the Kelvin-H elmholtz 
timescale of about 10 — 100 s (iPons et al.lll999D . and the 
wind should then become relativistic and non radiatively 
dissipative. A few days after the supernova explosion, the 
temperature of the soft thermal photons from the sur- 
face of the neutron star drops to T < 10^ K and photo- 
disintegration on this background radia tion can also be 
neglected, even inside the light cyli nder ([Protheroe et al.l 
[T9M IBednarek fc Prothel^[2nnl . 

In the pulsar wind beyond the light cylinder, pos- 
sible acceleration sites thus lie close to the equato- 
rial plane of the star, at a distance i?a > ^mXn = 
max(rmi„,iin,rmin,c,rmin,c,7r) ~ i x lO' ''^ Cm, assum- 
ing 7l ^ 10'^. The fact that r,„i,i > rdiss implies that the 
unipolar induction toy-model could apply, and that par- 
ticles could reach ultrahigh energies within this range of 
distances. 



2.3. Heavy nuclei injection 

One can mention three channels via which heavy ions 
could be seeded in the neutron star wind. Note that sce- 
narios of pulsar winds loaded with heavy nuclei give a 
satisfactory explanation to some observations. For in- 
stance, the morphological features of the Crab Nebula 
could be the signature of resonant scatte ring of pairs of 
electrons and positrons by h eavy nuclei (jHoshino et al.l 
[T99I IGallant fc Aronsl[T99l . 

The classical argument that applies best in 
our scenario is that iron nuclei can be stripped 
off the neut ron star surface, as has bee n sug- 
gcsted bv iRuderinan fc Sutherland! (|1975| ) and 
lArons k, ScharlemannI (|1979( ). Strong electric fields 
combined with bombardment by particles can extract 
ions from the polar cap regions, where the co-rotation 
charge is positive provided that • B < 0. The surface 
of a neutron star being composed mainly of iron-peaked 



elements, it is possible that heavy nuclei get injected in 
the wind by these means. 

Heavy nuclei loading of the pulsar wind by mixing 
of the stellar material via Kelvin-Hel mholtz instabilities 
or oblique shocks was also proposed (jZhang et al.l 120031 : 
IWang et al.l I2008D . This mechanism requires however 
that a jet goes through the stellar core, a case that is 
not considered in the present study. Kelvin-Helmholtz 
instabilities might also occur at the interfac e between 
the wind nebula and the su pernova remnant (jJunlll998l : 
Ivan der Swaluw et al.ll200^ . but it is unlikely that the 
envelope in that region has a metallicity high enough to 
mix large amounts of heavy nuclei in the wind. 

The nucleosynthesis of heavy elements by r-process 
in the neutrino-driven wind at the very early phase of 
the proto-magnet a r form ation has also been discussed 
by iMetzger et all ()2011allb[ ). These authors find that 
the production rate of nuclei with A > 56 can be 
important during the first 1 to ~ a few 100 s, when 
the electron fraction Yg could be fairly low, the wind 
expansion time Texp ^ 10^ s, and the entropy S < 
100 fcb nucl eon~i, as is required for a successful r-process 
(see, e.g., I Hoffman et al.l Il997l ). Though these results 
are obtained for the case of a highly magnetized proto- 
magnetar driving a jet (as in iBucciantini et al.l l2007| ) , 
they can be applied in a non-coUimated mildly magne- 
tized wind case, as the evolution of S and Toxp is mostly 
ruled by thermal ingredients (and the rotation speed) in 
the times considered. However, we will see in the next 
section that the supernova envelope at t ~ 10 — 100 s 
is too dense to allow the escape of particles, whatever 
their mass number. At later times, as the wind cools 
and becomes relativistic, the neutrino heating efficiency 
drops, shutting off the r-proccss. It is thus unlikely that 
this channel can seed heavy nuclei in the wind in our 
framework. 

3. UHECR ESCAPE FROM SUPERNOVA ENVELOPES 

Particles accelerated in the pulsar wind further need 
to escape from the pulsar wind nebula itself, and then 
from the surrounding young supernova envelope. We as- 
sume in this study that the supernova envelope is not 
totally disrupted by the wind, and that particles do not 
escape through a region punctured by a jet, like in a 
str ongly magnetized prot o-magnetar scenario discussed 
by IMetzger et al.l (|2011a[ ) — see Appendix |X] for further 
discussions. 

The escape of acceler ated ions from the magnetar wind 
nebula was discussed bv lArond (|2003f ). In Section [Z!2l we 
argued that at distances r ^ i?L , the curvature radius of 
the ions reads: pc ~ 2r]r^/Ri^ ^ r. Hence, particles are 
not coupled to the magnetic field lines and can escape 
the wind beyond Vnun- 

In supernova envelopes, m agnetic fields are of order 
a few mG at most (see, e.g., iRevnolds et al.ll2011l for a 
review). The Larmor radius of the ions is thus much 
larger than the size of the envelope and their trajectories 
can be treated rectilinearly. We give in the following 
section, estimates of the density profile and composition 
of young supernova envelopes, that we use to study the 
escape of UHECRs analytically and numerically. 



5 



3.1. Supernova envelopes 



As discussed for instance bv lChevalied ()2005f ). rotation- 
powered pulsars can originate in various types of core- 
collapse supernova: in Type II supernova resulting from 
red supergiant stars with most of their hydrogen enve- 
lope intact (SNIIP), or with most of their hydrogen lost 
(SNIIL and lib), or in Type lb or Type Ic supernova 
(SNIb/c) that stem from stars wi t h all their hyd rogen 
lost. See also iMaeda et al.l 120071 IWooslevI [2010l and 
IPiro fc Ottil201lL IK asen fc Bildstenll2oTor for supernova; 
associated with magnetars. iChevalieii (120051 ) finds that, 
of the remnants with central rotation-powered pulsars, 
the pulsar properties do not appear to be related to the 
supernova category. 

Within a few days after the explosion, the supernova 
enters a free expansion phase with velocity distribution 
V = r/t, that lasts several hundreds of years. A straight- 
forward way to model the evolution of the density of the 
ejecta is to assume that the ejected mass Moj will expand 
spherically in time with a mean velocity v^i over a shell 
of radius i?sN = Wej*- The ejected velocity, i?ej, relates 
to the supernova explosion energy and the ejected mass 
through: 
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where we defined Mej^io = A/cj/lOA/© and i?cj,52 — 
i?ej/10^^ ergs. Most core-collapse supernovas are inferred 
to have explosion energy i?ej ^ 10^^ ergs. However, for 
the pulsars with millisecond to sub-millisecond periods 
considered here, one can expect that the rotation en- 
ergy of order (l/2)/r2^ ~ 10^^ ergs will be transfered 
with in a fraction of a year to the surrounding ejecta 
(see iKasen k, BildstenI I2O1O0 . Depending on the radia- 
tion conversion efficiency of this energy, the surrounding 
supernova could become ultraluminous. Some ultralumi- 
nous SNIb/c and SNII have indeed b een detected with an 
explosion energy > 10^^ ergs (e. g .. iNomoto et al.l [20011: 



^, jrgs (e. g ., L 

WooslevI [2OIOI: I'Piro fc OttI l201lt iBarkov fc Komissarovl 
20100 . 



The mean density over Rs^{t) can then be written: 



/Osn(0 = 



2 X lQ-''Mlll,E;.'i^t-^ g cm-3 , 

(10) 

where tyr = t/\yi, which is the timescale to reach a 
pulsar spin that enables the acceleration of iron up to 
~ lO^O-S eV (see Eq. U). The column density integrated 
over i?gN as a function of time reads 



2/sn(0 — PsnRsn 
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More detailed modelings show that the density evo- 
lution of the ejecta is expected to depend on the type 
of supernova. Yet, we demonstrate in what follows that 
Eq. ([TT]) above provides a good estimate for the evolu- 
tion of the integrated column density of various types of 
supernova envelopes. Indeed, we will see in the next sec- 
tion that the escape of UHECRs is determined by their 
interactions on the baryonic envelopes. Because these 
interactions solely depend on the integrated column den- 



sity of the envelope, the detailed density profile is not 
crucial to our calculations. 

Under the assumption of adiabatic, spherically 
symmetric flows, the numerical calculations of 
iMatzner fc McKed ()1999D show that the density of 
a Type II supernova in the dense central region can take 
values as high as: 



PSNIl(t) 10 ^^A'4^/?o^cj,52^V^ 



gem 



-3 



(12) 



Most type II supernovae eject a mass of order Moj^io 
(jWooslev fc Weaveilll995D . This dense, relatively flat re- 
gion extends to radius i?b ^ 2(i?cj/Mcj)^/^t and is sur- 
rounded by a steep outer power-law profile. The column 
density that the accelerated particles have to traverse to 
escape is then: 

2/SNIl(t) = PSNII-Rb ^ 4A/^j_io£'^yg2*yr^ gcm"^. (13) 

For Type Ib/c/bc superno vae, one can apply the model 
of IMatzner fc McKed ()1999D for the explosion of a star 
with a radiative envelope, which yields: 



PsNib/c(i) = 7 X 10 



-17 



O.Olc 



-1.06 



p-0.97.-3 „„^-3 

^cj,52 V gem 



(14) 



out to radius i?b, beyond which the density decreases 
steeply. We have assumed in this estimate an explosion 
energy of £^ej,52 and an ejecta mass of Moj_2 = Afoj/2 Mq, 
which are derived from the observation of such objects 
(jProut et al.ll2010[ ). The corresponding column density, 
taking into account the velocity distribution v = r/t, 
reads 



Vsmh/cit) — I PsNib/. 
Jo 



Ar ^ 9 Ml,E-.%t;^' gcm-2 



(15) 

Equations pT|) . ([TS]). and ((T5|) agree within factors of 
a few. It is thus reasonable to consider Eqs. PH)) and 
pip as representative of the envelope mean density and 
column density, for types II and Ib/c supernova. Equa- 
tions and ((T5|) show that higher ejecta energy E^j 
and lower masses Mej would enhance the column den- 
sity. The effects of such cases on particle escape are also 
discussed throughout the paper. 

One can further note that if the pulsar wind shr eds its 
surrou nding supernova envelope, as discussed in lAronsI 
l)2003f ) for the magnetar case, disrupted fragments would 
expand in the interstellar medium. In this case, one can 
weight the initial supernova density by C~^^^, C = Sp/p 
being a factor meas uring the dumpiness of the envelope 
(|Murase et al.ll2009( ). A high C would ease the escape of 
UHECRs from the envelope. However, the values of C 
remain difficult to evaluate, as no observational evidence 
of such phenomena has been detected. 

The composition of the supernova ejecta depends upon 
the type, progenitor mass, and the final interior mass 
of the supernova. CXO J164710.2-455216's associa- 
tion with the Westerlund 1 star cluster argues that at 
least some pulsars arise from massive star progenitors 
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uno et ani2006D . But, as mentioned before, rotation- 
powered pulsars and magnetars have been invoked for 
a wide variety of supernova types. The composition of 
a type lb sup ernova is roughly 50% helium and 50% 
C/0: e.g., the IWooslevI (|2010l ) progenitor is ~50% he- 
lium, ~43% carbon and ~7% oxygen. Type Ic super- 
novae (more numerous than lb supernovae) are com- 
posed almost e ntirely of C/0 and heavier elements: e.g., 
iMazzah et all ()2010( ) argued that SN 2007gr was com- 
posed of roughly 75% C, 15% O, 8% ^GNi, and 2% S. 
Type II supernovae have a range of ejecta, ranging from 
roughly 60% H, 30% He, and 10% C/0 to explosions very 
similar to type lb supernova with small amounts of H. 

We will discuss in Section 3.4 how the escaped UHECR 
spectrum varies between pure hydrogen and pure helium 
envelopes (or helium and carbon envelopes). 



3.2. Analytical estimates 

In accord with the discussion at the beginning of Sec- 
tion [3?T1 we will consider in the following that Eqs. (|10|) 
and (fTTj) provide a reasonable estimate of the evolution 
of the density of the supernova envelope surrounding the 
neutron star. 

Successful escape of UHECRs from the envelope will 
occur if the shell crossing time tdyn is shorter than the 
cooling time by hadronic, thad, and photo-hadronic, t^^j, 
interactions. 



The acceleration of a particle to the energy E happens 
at a time after pulsar birth: t{E) ~ tspin{E). We can as- 
sume that the thickness of the supernova shell to traverse 
at a time t is given by i?sN — tvcj . Indeed, from the val- 



ues given in Section [521 the acceleration site i?a Rsn, 
as soon as t > 100 s. The crossing time for UHECRs 
traveling at the speed of light then reads: 



tdyn{E) 



-^tspin{E) 



(16) 



As the expansion time scale of the envelope is tc 



Rsn/Vc 



t. 



< t 



dyn: 



one can neglect the evolution 



of the envelope density during the escape of a particle. 

The timcscalc for hadronic interaction losses can be 
expressed as: 

th,d{E) = mi,{cpsN[tspin{E)]<j{E)^{E)}-^ , (17) 

where is the mass of the dominant target ion com- 
posing the envelope. The parameters £,{E) and cr{E) 
are the elasticity and the cross-section of the interaction 
at energy E. For our analytical estimates, we evalu- 
ate their valu es roughly from th e hadronic interaction 
model EPOS (j Werner et al.ll2006[) . We assume that the 
cross-sections of the hadronic interactions do not vary 
strongly above E > IQ-'^^ eV and set them to dp = 130 mb 
for proton-proton interactions and (Tpc = 1.25 b for iron- 
proton interactions. The number of nucleons carried out 
at each interaction can vary from 1 to A—1, with a large 
spread in values. For demonstrative purposes, we take 
an average value of ^ = 0.4 for both p — p and p—Fe in- 
teractions. These calculations are done accurately using 
EPOS in our numerical calculations in the next section. 



The condition of escape from the supernova envelope 
can thus be written as idyn/^had < 1, yielding 

t > icscp - 1.2 X 10^ A/cj,io£^ej,52^ « for proton, (18) 
t > icscFc - 3.8 X 10^ A/oj,ioSej,52' « iron, (19) 

where we assumed a supernova density profile following 
Eq. (flO)) . Cosmic rays at ultrahigh energy will escape 
only if they are produced at late times t ^ 1 yr, when 
the envelope density has decreased. Because nuclei of 
charge Z at a given energy E are produced at a time 
ispin(£^) oc Z (Eq. H]), one has tdyn/^had oc Z^^. The 
escape condition from the baryonic envelopes at a fixed 
E should consequently be eased for heavier nuclei. 

Still assuming the supernova density profile of Eq. (|10p , 
and using the spin-down time given in Eq. one can 
express the cut-off energy above which injected primary 
particles should not be able to escape the envelope: 



E, 



cut,Z 



y47rmbUej / 



9 /c3 



(20) 



1/2 



^ 7.5 X 10i« Z.iuh^ti^L^I^wE'J^h (^) ■ cV(21) 

^ 1.2 X 102"Z26rn/45M3~a5^^cTao£^S (v)'^'"^^^ 

where the first numerical application corresponds to pro- 
tons and the second to iron nuclei. Note that under the 
crude approximation that a oc A^/'^, £^cut,z oc Z/A'^/'^. 
For Ecut,z ^ Ei, Ecut.z does not depend on 

This trend is illustrated in Figure [U where the main 
timescales at play arc displayed: idyn and ihad as a func- 
tion of particle energy E, for various pulsar parameters 
r2 and and for both pure iron and pure proton injec- 
tions. As expected, iron particles can escape the enve- 
lope at higher energies, as they can reach these energies 
at later times. Lower magnetic fields (// < lO'^^) lead to 
longer tgpin at a fixed E (Eq. |4|) , while high pulsations 
{Vl > 10^ s) lead to higher acceleration energies (Eq. [1]). 

When iron nuclei are injected, secondary particles are 
produced by hadronic interactions for times t < tcsc,Fc- 
These secondaries of mass and charge numbers [A, Z) can 
escape the envelope only at times t > tosc,z, where tcsc,z 
is defined as the time at which idyn/^had = 1- Hence, 
secondaries that will escape from the envelope have nec- 
essarily been produced between iosc.z < t < tcsc,Fo, i-C, 
the lightest secondaries will escape first. This translates 
in terms of the energy range of the primary iron to: 
i?cut.Fc < E < Epc(tesc,z), where we can further express 
Epcitcscz) = (26/Z)i?cut,z- The main fragment among 
secondary particles will thus emerge from the envelope 
between energies E[^ 

A 



low.z < < E'z, with 



E\ow,Z = ^£'cut,Fo 



^2.1 X 10i«A,7iJ45M3"a5^^cTAo<52eV , 



^-'■^=56Z ' 



cut.Z 



3.5 X 10'^Ar„h,^,,,\,M-^^,,E:!;,eY 



(23) 



(24) 
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Fig. 1. — Timescales at play for the escape of UHECRs from a supernova envelope vfith A?fcj,io and ^'cj,52- The crossing time tdyn 
(dashed lines) and energy loss time by hadronic interactions thad (solid lines) are displayed as a function of particle energy E, for pure iron 



(red) and pure proton (blue) injections. The timescales are calculated for various pulsar initial rotation velocities Qi 
magnetic dipole moments /i = 10'^'' to lO'^^'^ cgs, as labeled. The other pulsar parameters are set to 7 = 10^^ gcm^, r) - 



= 10^, 

: 0.1. 



lO^s-i and 



The numerical estimates are calculated for secondary 
protons. Peaks of the various secondary elements should 
appear in the escaped cosmic-ray spectrum at their re- 
spective energies. A tail due to lower energy secondary 
nuclcons {E < E[^^ ^) following approximately a power- 
law with index ^ — 1/2 should also be produced together 
with the main fragment, down to PeV energies. The am- 
plitude of this tail around ^ E[^^ ^ is about a fraction 
of the number of the main fragment. 

From Eq. dH]) and (fTT]) . one can derive: 
tdyn/^had ^ 3 X 10^0 ^^^A'/^. i^i^.^^^gj at 1-2 = i/100 s, 
for {A, Z) = (90,40). We assumed a cross-section 
(jgo = 1.5 b for nuclei- proton interactions and an elas- 
ticity of ^ = 0.4, at energies E ^ 10^" eV (in the target 
rest-mass frame). This demonstrates that nuclei with 
^ > 56 that could be injected at times t ^ 10 — 100 s 
if a successful r-process occurred in the neutrino-driven 
wind (Section 12. 3p , cannot survive the crossing of the 
supernova envelope. 

Ultrahigh energy ions could also experience photo- 
disintegration in the radiation fields generated at the in- 
terface between the pulsar wind and the supernova shell. 

This radiation field can be expected to be significant 
if the supernova explosion is driven by the pulsar wind, 
as expected for millisecond rotators. A fraction 77^ of the 
wind energy ^ {1/2)IU? can be converted to radiative 
energy via internal shocks and another fraction 77th of 
this radiation then thcrmalizcs depending on the opacity 
of the medium. This thermal component peaks at energy 



time by photo-disintegration of order: 



tA'r,th = [cCAji^<^A-y/<^Aj)(^A-yUth/t 



'10^ A 



-0.21 

56 



£'cj,52 



3/8 



M 



-9/8,9/4 

oj.lO Vi- 



(25) 
(26) 



kT 



ergy density Uth ^ 
where rj^^i = yy^/O.l 



-3/4 



eV, 

3, 



with en- 

-3 



'oj,52 ^"oj.io'-y 

erg cm 

This background leads to a cooling 



where AeA-y/s A ^ ~ O AA^;?^, a a-, ^ 8 x 10"^*^ ^456 cm'^ 
(jMurase et al.l l2008| ). and we take for the elastic- 
ity of the ^7 interaction: S^A-y — ^/A (which 
is a crude approximation). This estimate of 
the cooling time is valid for cosmic-ray energy 

EA,pe.^ 4 X 1017 {Vf,mb)-'^^El.%M;j%^4i\v, and 
is about one order of magnitude larger for Ea ^ EA,pcak, 
as the photo-disintegration cross-section lowers. At the 
highest energies {Ea ^ 10^° eV), photo-disintegration 
could thus play a role on the escape of cosmic rays if the 
radiation and thermalization efficiencies are higher than 
ViVth ^ 10~^. The rate of wind energy going to radiation 
is eva luated to be of order 10% (e.g., iKasen fc BildstenI 
l201Clf l. but the thermalization fraction of these photons, 
77th J is not known, due to the uncertainties on the 
opacities in the internal shock region. Mixing and 
Rayleigh- Taylor instabilities effects creating finger-type 
structures could lead to a leaking of the high energy 
photons, and the thermalization fraction could be as 
low as < 10%. A higher acceleration efficiency 77 would 
also enable particles to reach the highest energies by the 
time the radiation field intensity has become negligible. 
Given these uncertainties, and for simplicity, we will 
assume in this paper that the radiation field can be 
neglected for the escape of UHECRs from supernova en- 
velopes, the baryonic background playing the major role. 



To summarize, the conditions for successful accelera- 
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tion and escape above 10 eV can be written as: 



Bnf > (1012-4 G) X {10^ s-^f Z^qIj-^R 



i 1 /2 
'7l-'45^'47,10^oj,52-^*,10 



^1 D-3 
-3 



(27) 



Higher values of the magnetic field would allow higher ac- 
celeration energies, but would require lower ejecta mass 
and higher explosion energies. Note that 10 Mq can 
be viewed as an up per bound for the e jecta mass for 
type II supcrnovse (jWooslev et alJ l2002| ) . One might 
also advocate that the presence of clumps could lower 
the overall densities and allow the escape of particles at 
E > lO^" eV. All in all, the parameter space allowed 
for successful acceleration and escape appears to be nar- 
row, but we will see in Section 14.11 that the low rate of 
sources required to account for the observed UHECR flux 
would compensate for this issue. A higher acceleration 
efficiency rj would also broaden the allowed parameter 
space. 



3.3. Numerical Setup 

As discussed in the previous section, the hadronic in- 
teraction between UHECRs and the baryonic envelopes 
is the determinant factor that would affect the injected 
UHECR spectrum. 

The interactions with the baryonic envelopes were cal- 
culated by Monte- Carlo for injec t ed nu clei and their 
secondaries. As in iKotera et al.l (I2009D. we used the 
hadronic interaction model E POS (iWerner "eFall [20061) 
and the fragmentation model of lCampi &: Hiifned ( 198l[ r. 
as imple mented in the air s hower simulation code 
CONEX (|Bergmann et al.l[2007[) . 



In the case of a non-hydrogen baryonic envelope, the 
interaction products can be derived from the nuclei- 
proton interaction case by a superposition law. In the 
target rest frame, the products of the interaction between 
a projectile of mass number and energy ( Aproj , ^'proj ) and 
a target nucleus of mass number ^targ are roughly equiv- 
alent to Atarg times the products of the interaction be- 
tween a projectile with (Api-oj, i'proj/^targ) and a target 
proton. The exact cross-sections are nonetheless com- 
puted with EPOS. 



In the simulations, 
angular velocity fii ^ 
ment /i ~ 10^°"^^ cgs, 
netic dipole field B 
is an upper limit (^ 



we modeled pulsars with initial 



10 



3.0-4.2 



s and magnetic mo- 



corresponding to a surface mag- 
2 X 1012-15 G. Notice that there 

10^-2 s~i) on the initial angular 

speed (jHaensel et al.l |1999() . For each set of parame- 
ters, 10^ cosmic rays are injected following a power-law 
energy spectrum as in Eq. ([2]) with minimum injection 
energy E^in = 10" eV, and the maximum acceleration 
energy Ei calculated in Eq. Above Ei, the spec- 

trum cuts-ofF exponentially. Nuclei with initial energy 
E are injected at a radius i?a = IQi*' cm (correspond- 



10'* 



see discussion in See- 



ing to ^ 3 T„un for Hi 
tion 12. 2p at the time tspin(£'), and propagate through a 
supernova envelope of total ejected mass 10 Mq (2 Mq 
in Type Ib/c supernova case) expanding at a constant 

rate Vej = 10^ ii^Vg^MT^p^ cms~i. The evolution of the 
ejecta density is assumed to follow Eq. (fTO|) . We studied 
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Fig. 2. — UHECR spectrum before (dash) and after (solid) escape 
from a hydrogen supernova envelope with Mej lo and -Eej,52, with 



pure proton injection. The pulsar pa 
rf = 0.1, n = 10*-0 s-1, and fi = 10^° 



arameters are / = 10 gem , 

-5 ncya 



pulsars embedded in pure hydrogen, helium and carbon 
supernovas. 

3.4. Numerical Results 

We first assume a pure hydrogen envelope. The re- 
sults are presented in Section 13.4.11 Simulations using 
more supernova envelopes with heavier composition are 
discussed in Section [3.4.21 

3.4.1. Pure Hydrogen supernova envelope 



Figure [2] presents the injected (in dash line) and es- 
caped (in solid line) spectra of pure proton injection by 
a pulsar with initial angular speed 51i = lO'^s-i and 
magnetic dipole moment fi = lO'^'^'^cgs. The injected 
spectrum follows the characteristic —1 spectral index in 
Eq. ([2|). As predicted in Eq. ([2T|) UHE protons above 
^ 10 EeV fail to escape the supernova envelope, since the 
region is still very dense at the time they are produced. 
Below a few EeVs protons are free to escape. Protons 
with energy in between can partially escape with signifi- 
cant flux suppression. EPOS shows that for one 10 EeV 
primary proton, the peak of interaction products lies at 
IQi^ eV; the chance of resulting a secondary proton with 
E > lOi^eV is less than 0.01. Therefore we can barely 
see the secondary protons in our energy window of sim- 
ulation. 

The spectra of pure iron injection by pulsars with 
Qi = 10^ s-i and = lO^O-^ lO^i-^ cgs are shown in Figd 
In the top plot (^ = 10^°-^ cgs, ft = lO^s"!), primary 
iron nuclei with energy up to i?cut,Fc = 1-2 x 102° eV 
can escape without significant loss. As discussed in our 
analytical estimates, most secondaries should originate 
from primary iron nuclei with energy between i^cut.Fc ~ 
1.2 X 1020 eV and 56 x E^ut.p = 4.2 x 102° eV, correspond- 
ing to the iron cutoff and iron mass number times the 
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Fig. 3. — UHECR spectrum before (dash) and after (solid and 
dash dotted) escape from hydrogen supernova envelope v^ith M^j ig 
and -Eoj,52i with pure iron injection. The pulsar parameters are 
I = 10*'5gcm2, ri = 0.1, Q = 10'' s^^, and fj. = lO^O '^cgs (top), 
and fi = lO^'^-^ cgs (bottom). Different compositions are listed as 
in the legend box. 

cutoff of secondary protons. In agreement with Eq. (j24p . 
secondaries lie between (1.0 — 5.0) x 10^^ eV for proton, 
2.0 X lO^s - 1.3 X 1019 eV for helium, 7.9 x lO^* - 4.0 x 
10^9 eV for CNO, (1.3 - 7.1) x lO^^eV for Mg-like ele- 
ments and 2.0 x 10^^ - 1.1 x 10^° eV for Si-like elements, 
with the peak positions scaled to the mass number of the 
elements and the bump width being almost the same in 
logarithmic coordinates. The significant tail of protons 
below 1 EeV comes from the products of the hadronic 
interactions. On average, each interaction of a 500 EeV 
iron nucleus results in one EeV proton among its prod- 
ucts. The strong signals from secondary nuclei contribute 
to a steeper overall spectrum (in solid black) which fol- 
lows ~ E-'^ at 10^^-^ - 102"cV. 



When the magnetic field is 10 times stronger (/i = 
lO^i-S cgs, 17i = lO^'s-i, bottom plot of Fig [3]), the pul- 
sar spins faster and the cutoff for primary and secon- 
daries are lowered by 10 times (sec Eq. [M]) . Hence, 
the ij. = lO'^i'^cgs case presents a similar shape as the 
/i = 10^°-^ cgs case except an overall shift to lower ener- 
gies by a factor of 10. 

As pointed out in Section 12.11 at low energies when 
E ^ Eg the gravitational wave losses are negligible and 
tgpin is independent on the initial rotation speed fii for 
E Ei. A pulsar with higher initial angular veloc- 
ity can inject UHECRs with greater maximum energy. 
However a mi nimum spin period ^ 0.4 ms is allowed for 
neutron stars (jHaensel et al.l[T"999( ) corresponding to an 
upper limit (^ 10*'^ s"^) on the initial angular speed. 
Magnetic dipole moments iJ. greater than 10'^^ cgs would 
make the spin-down process too fast to allow UHECR 
escape. On the other hand pulsars with ^ < 10"^" cgs are 
not energetic enough to accelerate particles to ultrahigh 
energy (see Eq. [27]) . To determine the best escaping re- 
gion we ran a parameter scan with 15 x 15 sets of (17, fi) 
and the results are presented in Figure HI 

We define the cut-off energy Ecut as the energy the ra- 
tio between the escaped and injected particles is less than 
10%. It corresponds approximately to the highest energy 
of escaped cosmic rays -Ecut.z defined in Eq. (|2T|) . In Fig- 
ure 131 the contours represent i?cut reached after escaping 
hydrogen supernova envelopes with Mej^io and i5ej,52 for 
pulsars with dipole moment /i and initial angular velocity 
fli. In the proton case (top), protons with energy above 
10^0 eV cannot escape the supernova envelope in our 
model. In the iron contours (bottom) , the parameter re- 

gion with ifl « 1030-00-30-72 cgs) X {fti « 103-95-4.20 

allows the escape of iron nuclei with energy greater than 
10^° eV. This parameter scan is based on a supernova en- 
velope with density profile described in Eq. ([T0| . Higher 
values of explosion energy and lower ejecta mass could 
lead to a broader enclosed parameter region that allows 

the escape, as E'cut.z scales with M^^\qEI-'^2 (Eq- BUI) . 
Our res ults agree with t he theoretical prediction from 
Fig.l in lBlasi et aD(|2000( ). except that we have a smaller 
parameter area that allows escape. This comes from our 
assumption that only rj ^ 10% of the induced potential 
turns into UHECR energy. 

3.4.2. Helium-Carbon/Hydrogen-Helium supernova 
envelopes 

Results in Figure [5] are from simulations with an ejecta 
mass of Moj = 10 Mq, explosion energy E^j = 10^^ ergs, 
and a composition of pure ^He (top) and pure ^^C (bot- 
tom). As discussed in Sec. 13.31 realistic envelopes for 
SNII and SNIb/c are more complicated and could be 
evaluated by a combination of Fig. [3] and [H Spectra of 
UHECRs escaped from envelopes abundant in heavier 
elements maintain features from that with a pure hydro- 
gen envelope. For instance, in case of a pure helium 
envelope (top plot in Fig. [5]), the spectrum preserves 
the 'original' secondary peaks at 6.3 x lO^^eV for hy- 
drogen, 2.5 X 10^^ eV for helium, 8.0 x lO^^ cV for CNO, 
1.3 x 10^9 cV for Mg-like elements and 2.5 x 10^9 gy for 
Si-like elements. These peaks are similar to the ones in 
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Fig. 4. — Parameter space with cut-off energy (Ecut.z) contours, 
for a hydrogen supernova envelope with Afoj.io and ^'ej,52i a^nd 
pulsar parameters 745 and ri\ . The solid lines refer to cut-off par- 
ticle energies after the escape. Up is proton injection and down is 
iron injection. Notice that current neutron star models suggest an 
upper limit of rotational speed at f2i < lO^ '^s^^. Note also that 



^'cut z scales with M 



-1 



^1/2 
cj,10-^cj,52 



(Eq.[20ll. 



a pure hydrogen envelope (see the first plot of Fig. 
except that they are located at 4 times lower in energy, 
due to the 4 times heavier interactant. 

The case of heavy envelopes can generate multiple 
peaks to the left of the original peaks due to multiple 
products. According to the superposition law, the num- 
ber of products scales with after N interactions with 
envelope baryons of mass number A. So the later gener- 
ations (tcrtiaries and so forth) whose energy are mostly 
below 10^* eV are far more numerous than the earlier 
generations (primaries and secondaries). This brings the 
low end of the original peaks up to be a second, or even 
third additional peaks for all compositions; they also con- 
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Fig. 5. — UHECR spectrum after escape from a supernova enve- 
lope with Mgj 10 and Eci,52, with composition (top): 100% "^He and 
{bottom): 100% ^^C. The pulsar parameters are / = lO^^gcm^, 
Tf = 0.1, Q. = 10* s~\ and fi = cgs. 

tribute to an increment of primaries around 4 x 10^^ cV 
for helium envelopes and 1.3 x 10^^ eV for carbon en- 
velopes. 

4. IMPLICATIONS FOR THE SCENARIO OF UHECR 
PRODUCTION IN NEWLY-BORN PULSARS 

The success of a UHECR source scenario lies in its abil- 
ity to reproduce these observations: i) the energy spec- 
trum, ii) the composition, iii) the anisotropy, and iv) on 
the fact that it requires a rate of sources consistent with 
the population studies inferred from other astronomical 
observations. 

As wc discuss in this section, the results obtained in 
this paper suggest that all four points could be rea- 
sonably achieved in the extragalactic rotation-powered 
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pulsar scenario. Newly-born pulsars are natural candi- 
dates to reproduce points ii) and iii), due to their iron- 
peaked surface (if the composition at the highest ener- 
gies proves to be actually heavy, as the measurements 
of Auger seem to indicate) and their transient nature. 
Though point i) is challenged by the fact that the toy 
model of unipolar induction generates a hard spectrum 
that does not fit the observed UHECR spectrum, our 
results show that the slope could be naturally softened 
during the escape from the supernova envelope (also seen 
in Bcdnarck & Bartosik 20041 for Galactic pulsars). The 
range of parameters for the pulsar and its surrounding 
supernova allowed for a successful acceleration and es- 
cape at the highest energies is relatively narrow. This 
potential issue is however compensated by two advan- 
tages. First, the range of values required for the initial 
parameters of both the pulsar and its supernova are close 
to the ones inferred for the youngest is olated pulsars ob - 
served nowadays (see e.g., Table 3 of iChevalierl 120051 ). 
Second, the rate of such objects required to account for 
the observed flux of UHECRs is low, of order fs < 0.01% 
of the 'normal' (as opposed to binary millisecond) pul- 
sar birth rate. Point iv) can hence also be deemed as 
reasonably satisfied. 

We will also examine in what follows, the implications 
of our results on the arrival directions of UHECRs in the 
sky, and on possible probes of this source scenario. We 
also discuss the signatures expected for secondary mes- 
sengers such as neutrinos, gamma-rays and gravitational 
waves. 

These implications are first discussed under the as- 
sumption that the currently observed UHECR flux has 
an extragalactic origin. The contribution of Galactic pul- 
sar births is discussed in Section [ 



4.1. Required source density and type of source 

The magnctar birth rate necessary to a ccount f or the 
observed flux of UHECRs was esti mated inlAronsI (|2003f ) 
and updated for various cases bv lKoteral ([20Tl[ ). The 
same calculations can be applied to our case, for rotation- 
powered pulsars. 

For a population of identical neutron stars with ini- 
tial rotation velocity fij and magnetic dipole momentum 
/i, satisf ying Eq. (l27l) . one can adapt the normalization 
found by iKoteral ()2011h for negligible gravitational wave 
losses. An identical neutron star assumption is accept- 
able, in so far as the allowed parameter range of sources 
for particle acceleration and escape is fairly narrow 
(Eq. [221). A birth rate of n ~ 10"* fisiZ^^} Mpc'^ yr'^ 
is required to produce the observed UHECR flux, in the 
absence of source evolution history. When the emissivity 
of UHECR sources is assumed to follow the star forma- 
tion history, the pulsar birth rate at z = is of order: 
nsFR 0.8n ~ 0.8 x 10" Vai Mpc'^ yr"!. This cal- 
culation assumes that the total Goldreich- Julian charge 
density is tapped in the wind for UHECR acceleration 
(Eq. [2]) . A lower efficiency would result in a lower energy 
flux per source, and thus in higher required densities. 

The above rates correspond to a fraction /g < 0.01% 
of the birth rate of 'normal' pulsars, which is of or- 
der 1.6 X 10~^ Mpc""^ yr~^ (or one per 60 years in 



the Galaxy, w hich is consistent with the supernova rate, 
lLorimeJl2008D . Among the 'normal' pulsar population, 
it is difficult to infer the number of objects that would 
satisfy Eq. (|27l) . as the distribution of pulsars according 
to their initial rotation velocities and magnetic field is 
not straightforward ( see examples of models discusse d in 
iGiller fc LiTOkH 12003 and iBednarek fc Bartosik|[200l . 

iFaucher-Giguere fc Kaspil (|2006f l find that the birth 
spin period distribution of pulsars is normal, centered 
at 300 ms and with standard deviation 150 ms, and that 
the initial magnetic field follows a log-normal distribu- 
tion with (log(B/G)) - 12.65 and triogs ~ 0.55. They 
stress however that this distribution of birth spin peri- 
ods is not precisely constrained by their method, and 
considerable deviations from this statistics could be ex- 
pected. Such a distribution would imply that < 2% of 
the 'normal' pulsar population could be endowed with 
sub-millisecond periods at birth. 



Equation ([27|) further depend on the supernova char- 
acteristics (ejected rnassand energy). However, as dis- 
cussed in iChevalieil (|2005l ) , the pulsar properties do not 
appear to be closely related to the supernova category. 
This introduces an additional degeneracy on the type and 
total number of objects meeting the requirements for ac- 
celeration and escape. Nevertheless, it is promising that 
the range of values required for the initial parameters of 
both the pulsar and its supernova are close to the ones 
inferred for th e you ngest observed isolated pulsars (see 



mierrea tor tn e you ngest oDser 
e.g.. Table 3 of IChevalier 20051) 



Hence, /g is a small enough fraction to leave reasonable 
room for poorer injection efficiencies, and to account for 
the narrowness of the parameter range of Eq. ()27|) . 

One should also keep in mind that both HiRes 
(lAbbasi et al. 2009') and the Pierre Auger Observatory 
(jAbraham et al.. ,2010bi) report systematic uncertainties 
of order 20% on the absolute energy scale of the spec- 
trum, which should be considered for the evaluation of 



The distribution inferred by iFaucher-Giguere fc Kaspil 
(|2006( ) implies that pulsars with birth periods ~ 300 ± 
150 ms are about e '~ 30 times more numerous than the 
submillisecond ones. Such pulsars could potentially ac- 
celerate iron up to E{P = 100 ms) ~ lO^*' eV (Eq.[l])'. For 
extragalactic pulsars with a similar acceleration mecha- 
nism to the case we discuss here (i.e., only a fraction fs 
of the existing population leading to cosmic-ray produc- 
tion), the amplitude of the injected spectrum at these 
lower energies is well below the observed one (even if 30 
times more numerous, the hard ^ power-law, 
below the peak due to the secondary protons, only over- 
takes the observed ~ spectrum closer to ankle ener- 
gies). At these low energies, the diffusion of cosmic-rays 
in the intergalactic magnetic fields would further prevent 
them from reaching us, if the sources are located at tens 
of megaparsec distances. On the other hand, a Galactic 
population of these more numerous slower pulsar births 
may give important contributions to the cosmic ray spec- 
trum below the ankle (see, e.g., [Bednarek fc Bartosi3 
[200l . 
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4.2. Propagated escaped energy spectrum 

The cosmic ray spectrum observed by the Pierre Auger 
Observatory can be described as a broken power-law, 
E~^, with spectral index A ~ 3.3 below the break (called 
"ankle") around lO^^-^ eV, and A ~ 2.6 above, followed 
by a fl ux suppression above ^ 10^^-^ eV (jAbraham et al.l 
l2ni0bD . 

O ne issue of th e model advanced bv lBlasi et al.l (|2OOO0 
and lAroiii (|2003[ ) for the acceleration of UHECRs in pul- 
sars and magnetars is the hardness of the produced spec- 
trum, that hardly fits the observations described above, 
even after propagation. These models were introduced 
in the "AGASA era" , to account fo r the absence of GZK 
cutoff in the observed spectrum (|Takeda et al.l I1998D . 
They aimed at producing a hard spectrum (of spectral 
index —1, see Eq. [3]) to fit the highest energy end of 
the spectrum, beyond E > 6 x 10^^ eV, and do not 
fit the slope at lower energies. The latest experiments 
report however that a suppression reminiscent of the 
GZK cut-off is pres ent at the high est e nergy end of the 
UHEC R spectrum (jAbbasi et al.l 120081 : lAljraham et~all 
l2010bD . Hence, a hard spectrum need no longer be advo- 
cated to explain the measurements, and now constitutes 
a disadvantage. 

iKoteral (|2011[ ) proposed to alleviate this issue by intro- 
ducing a distribution of initial pa rameters of magnetars 
among their population (see also iGiller fc Lipskil ( |2002D 
for the Galactic pulsars case). Such a distribution re- 
sults in a distribution of the maximum acceleration en- 
ergy, and adequate values can be found to soften the 
integrated spectrum and fit the observations. The same 
calculation can be applied to the case of rotation-powered 
pulsars. 

Note also that in order to have the monoenergetic- 
type acceleration spectrum given in Eq. ([2]), the wake- 
field acceleration which is based on the ponderomotive 
force requires the magnetic field to be coherent over the 
acceleration region. However, much smaller coherence 
scales can be naturally expected, leading to a stochastic 
acceleration, that could also produce a E~'^ spectrum. 
Suc h cases have been studied in di f ferent contexts by 
e.g.. lChen et"aI1 (|200l : IChang et al.l (|2009h . 

The results of iBednarek &: Barto"si3 ()2004[) already 
show that the injection of iron nuclei and their escape 
through the pulsar nebula can lead to a softer spectrum 
due to the production of secondary nuclei. This feature 
was however not deeply discussed and highlighted, as the 
AGASA energy spectrum available at that time highly 
differed fro m the current observati o ns. B esides, the cal- 
culations of IBednarek fc BartosikI (|2004[ ) are based on 
simplified hadronic interaction cross-sections, and on the 
raw assumption that one interaction leads to the frag- 
mentation of the primary nucleus in two nuclei with dif- 
ferent mass numbers. 

Our detailed analysis demonstrates that, within the 
range of pulsar and supernova envelope parameters given 
in Eq. (P7)) and Fig. |31 the injection of heavy nuclei and 
their escape from the envelope naturally enables the soft- 
ening of the energy spectrum to indices of order ^ 1.5 — 2 
(Figs. [3l [5|). As explained in Section \3A[ this softening 



stems from the abundant production of secondary nucle- 
ons, helium and intermediate nuclei at low energies. 

After propagation and interactions in the intergalac- 
tic medium, the injection of particles at the source with 
index ~ 2 is expected to provide a good fit to the ob- 
served UHECR spectrum. Our escaped composition can 
be identified with t he mixed c omposition in t roduc ed by 
lAUard et ail ()2008[ ) (see also lAloisio et all 120111 ) that 
contains 30% of iron and as sumes a max i mum proton 
energy of Ep^^^^ - lO^^eV. lAUard et al.l (pOOl ) calcu- 
lates that an injection index of order 2.0 — 2.1 is required 
to adjust the observed UHECR spectrum after propaga- 
tion through the intergalactic medium. If one assumes 
that the source emissivity in UHECRs has evolved ac- 
cording to the star formation r ate, the required inj ection 
index at the source is of ~ 1.2 (jKotera et al.ll2010l ). 

The bumps and irregularities apparent in the escaped 
spectra (Figs. [3l [5|) should be attenuated by the prop- 
agation, a possible distribution of neutron star charac- 
teristics (essentially a distribution of the dipole moment 
/i, initial spin Qi), and especially the envelope chemical 
composition. 

The flux of particles with energy below the ankle 
should not overwhelm other (possibly Galactic) compo- 
nents. Our calculations show indeed that the escaped 
spectrum should become harder below E ~ 10^^ eV, with 
a slope of order —1.5 due to the tail of secondary protons. 
The flux of these lower energy particles should also be di- 
luted by the large dispersion of their arrival times, after 
propagation in the intergalactic and Galactic magnetic 
fields. 

The injection of a pure proton composition by neutron 
stars is likely only viable in models where the envelope 
column density is thinner by many orders of magnitudes 
compared to classical supernova; at early times. In this 
situation, the resulting UHECR observable quantities 
are similar to what has been discussed until now: a 
hard spectrum injection should be expected after escape 
(spectral index —1), that could be reconciled with the 
observed spectrum by in voking a distrib ution of neutron 
star characteristics, as in IKoteral (|2nTl . 

One probe of this scenario (both in the proton or iron- 
rich injection cases) would be a sharp cut-off of the en- 
ergy spectrum at energies above £'cut,Fc (or E'cut.p for 
pure proton injection). A mild recovery is indeed ex- 
pected if the maximum acceleration energy were E > 
-|^q20.5 i-jj^g observed cut-off in the spectrum would 

then be due to the GZK effect. 

4.3. UHECR composition 

Recent measurements by the Pierre Auger Observa- 
tory indicate that the cosmic ray composition transi- 
tions from being dominated by protons below the an- 
kle (~ 10^*'^ eV) to being dominated by heavier nuclei 
with average masses s imilar to Si or Fe at ^ lO^^eV 
(jAbraham et al.|[2010al) . The instruments located in the 
Northern hemisphere, HiRes and Telescope Array, seem 
to observe a light composition up to the highest ener- 
gies, though the results of the former remain consistent 
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with those of Auger within errors. We caution further- 
more that the composition measured by these experi- 
ments concern energy bins below ^ 4 x 10"'^^ eV, due to 
the lack of statistics at the highest energy end. 

Neutron stars are one of the most likely places to inject 
heavy nuclei abundantly, as we discussed in Section 12.31 
It is interesting to notice that the escaped composition 
resulting from such an injection indicates a transition 
from light to heavy nuclei around the energy observed 
by Auger (Figs. El [5]). 

Note that iBednarek fc BartosikI (|2004D found a sim- 
ilar transition, but did not devote much discussion on 
that feature. That finding was not necessarily appealing 
during the AGASA era, when the composition was be- 
lieved to be light at the highest energies. To account 
for the con tinuation of the flux above GZK energies, 
IBednarek^ Bartosik (2004) added to their Galactic pul- 
sar population, a pure proton extragalactic component, 
which lightens their overall composition at the highest 
energies. 

As mentioned in the previous section, our escaped com- 
position is similar to the low ii^p.ma x mi xed composi- 
tion in troduced by lAllard et al.l (|2008[ ) and lAloisio et al.l 
(|2011| ). The resulting composition after propagation in 
the intergalactic medium when such a composition is in- 
jected is shown to conserve the tra nsition between lig ht 
to heavy elements around 10^^ eV (jAUard et al.ll2008l ). 

The injection of a mixed composition with < 10% of 
iron would remain consistent with such a transition. In- 
deed, Figs, m E] show that the rate of secondary protons 
is more than 10 times higher than the rate of injected 
iron. Injected protons would cut-off below i^cut.p and 
would not overwhelm the escaped iron flux at the high- 
est energies. 

One can note that, depending on the detailed transi- 
tion from extragalactic to Galactic component, the com- 
position found here may induce a n anisotropy signal at 
lower energies as was discussed in iLemoine fc WaxmanI 
(j2009( ). With such a dominant heavy composition at 
ultrahigh energies (> lO^'^^'^eV), one expects that any 
anisotropy signal at the highest energies would have a 
similar structure around 2 EeV where the composition 
is proton dominated, about two times stronger. Such 
an anisotropy is not observed by the Auger observatory 
(jAbreu et al.ll20lia| ). which may question the composi- 
tion of the mild anisotropy found at the highest ener- 
gies or imply a more complex composition structure both 
for the extragalactic as well as the Galactic component 
around EeV. 

The injection of a pure proton composition is not ruled 
out either in our scenario, but is only favored under strin- 
gent conditions on the early envelope density. 

4.4. Distribution of events in the sky 

The radio. X-ray and gamma-ray signals of rotation- 
powered pulsars and magnetars are too weak to allow 
their detection beyond our Local Group. For this reason, 
a direct spatial coincidence between a neutron star and 
UHECR arrival directions is not expected to be observed, 
if the source is not born inside our Local Group. 



Nevertheless, the distribution of UHECR events could 
follow the large scale structures, where neutron stars 
should be concentrated. In particular, these objects 
should be frequently found in star forming galaxies. Such 
distributions would be apparent only if the deflections 
experienced by particles in the Galactic and intergalac- 
tic magnetic fields are small. Moreover, anisotropic sig- 
natures would only be distinguishable for ensembles of 
particles with the highest energies. Above E ~ Eqzk = 
6 X 10^^ eV, the horizon that particles can travel with- 
out losing their energy is limited to a few hundreds of 
megaparsecs and the distribution of sources in that local 
Universe appears anisotropic. 

Neutron stars can be considered as transient UHECR 
sources. Cosmic rays with energy above Eqzk can 
indeed only be produced during the first Ats ~ 
4 ^26'7i-^45A^30^5 yr after the birth of the neutron star. 

This implies that if secondary messengers such as neu- 
trinos, gamma-rays, or gravitational waves were pro- 
duced at the same time as UHECRs, they would not be 
observed in temporal coincidence with the latter. The 
time delay experienced by UHECRs in the intergalactic 
magnetic field is indeed of order ~ 10^ yrs for one degree 
deflection over 100 Mpc, which is much longer than the 
duration of the UHECR production. 

Transient sources could lead to bursts of events in the 
sky if the number of event per source is important, and 
the arrival times of particles is not dil uted by the disper - 
sion induced by magnetic deflections (jKalli et al.ll20li[ ). 
For extremely high energy protons and low intergalatic 
fields, the total dispersion time due to magnetic deflec- 
tion, Et, can be shorter than the detector exposure time, 
Toxp, the number of sources contributing to the observed 
UHECR flux inside a radius of / = 100 Mpc is of order 
Ns ~ {4:Tr/3)l^hTcxp ~ 0.4, using the pulsar birth rate 
inferred in the previous section, and Tcxp = 10 yrs. The 
number of events that can be detected from each source 
is: TVev = fuHECR^oxp/(-BGZK 477^2) ^ 2 X 10^, where 
we assumed a detector exposure of Acxp = 3000 km^, 
as for the Pierre Auger Observatory, and the cosmic-ray 
energy output per source f uhecr ^ 5 x 10^^ erg in our 
milli-second pulsar scenario. It is likely however that 
cosmic rays arriving from most directions in the sky ex- 
perience a significant dispersion in their arrival time, due 
to magnetic fields: Et > Texp. This should be the case for 
iron nuclei, unless one assumes unrealistically low mag- 
netic fields. In that case, the number of events detected 
from one source would b e reduced by a factor Tny n/Et. 
The reader can r efer to iKotera fc Lemoin^ l|2008f l and 
iKaUi et all (|2011[ ) for detailed discussions on the depen- 
dence of Et on magnetic field parameters. 

A direct identification of the source could be possible 
if a pulsar was born inside our Galaxy, or close enough to 
allow X-ray or gamma-ray observations. The dispersion 
of arrival times inside our Galaxy tTcai reads: 
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Here we noted Bturb and Aturb the turbulent magnetic 
field intensity and coherence length respectively, and / 
the distance of the source. The time delay Jicai expe- 
rienced by particles due to the turbulent Galactic mag- 
netic field is typically much larger than (jQai, due to the 
additional deflection due to the regular magnetic field 
component. 

For proton injection and a weak regular magnetic field 
component, this implies that if such an event were to 
occur, a burst in UHECRs with a typical rise and decay 
timescale of a fraction of year would be observed in the 
sky, from a time <5<Gai after the onset of the explosion that 
triggered the birth of the fast-spinning neutron star. In 
this case, secondary messengers propagating rectilinearly 
would also arrive at a time StQ^i before UHECRs. 

For iron nuclei injection, the highest energy elements 
come out of the envelope as heavy nuclei. These 
should reach the Earth after a time delay of <5<Gai 
1750yrs//^Qi^pj,. For very close-by sources (e.g., at 2 
kpc), Jicai could be of order Texp, leading to a sudden in- 
crease in the detection of ultrahigh energy events (about 
A^ev ~ 4 X 10^^ over Jioai ^ 70 yrs, for a source at 2 kpc). 
In that case, the time order of escape of the different 
chemical elements from the envelope should be washed 
out by the fact that Jicai > Aig. 

Particles at energies E < Eqzk should arrive with 
more consequent time delays, so potentially from young 
rotation-powered pulsars that arc detected nowadays. 
The dispersion in time should however be as consequent, 
and such events should not be detected as bursts, but 
only as continuously arriving particles. No spatial clus- 
tering from such sources is expected either, as the de- 
flections experienced by particles at these low energies 
should be large. 

If EeV or higher energy neutrons were produced by 
these objects, by interactions of accelerated nuclei in the 
envelope for example, they would propagate rectilinearly 
and would appear as point sources. However, nearly no 
time delay between the detection of the birth of the neu- 
tron star and the arrival of the particles is expected. Spa- 
tial correlations between pulsar positions and neutron 
events are thus expected only if a new birth actually oc- 
curs in the Galaxy. 

4.5. Secondary messengers 

The propagation of UHECRs in the intergalactic 
medium should lead to the production of cosmogenic 
neutrinos and gamma-rays by interactions on the Cos- 
mic Microwave Background. The expected cosmogenic 
neutrino and gamma-ray fluxes depend mostly on pa- 
rameters inherent to cosmic-rays themselves (their com- 
position and overall flux), but also on the injection index 
at the source and the sou rce emissivity evolu tion history 
for diffuse fluxes (see e.g.. lKotera et al.ll2010l for a param- 
eter scan over these astrophysical variables) . The cosmo- 
genic gamma-ray signatures further depend on the struc- 
ture and strength of the intergalactic magnetic flelds, be- 
cause of the pair production/inverse Compton cascading 
of photons in the intergalactic medium. 

For a source evolution following the star formation 



rate, as can be expected for neutron stars, an injection 
of pure proton or proton-dominated compositions with 
power-law spectral index ^ 2.0 — 2.5 would successfully 
flt the observed UHECR spectrum. The resulting dif- 
fuse cosmogenic neutrino flux would lie within the gray 
shaded region of Fig. 9 of iKotera et al.l (j2010[ ). For an 
iron dominated injection up to a few times 10^" eV and 
a proton dominated injection below 10^^ eV (as we get 
in Figs. [HI E]) one expects a lower neutrino flirx, peak- 
ing around ~ lO^-^ GeV with E'^{dN/dE,,)\m^^ - 
5 X 10^ G eV cm~^ s~ ^ sr~^ (red dash-dotted line of 
Fig. 7 of IKotera et al.l l2010l ). For the diffuse cosmo- 
genic gamma-ray background, the same fit to the ob- 
served UHECR spectrum gives fluxes peaking around 
E^ ~ 10 GeV of order E?^{dN/dE^)\,n^^ - 7 x lO^^^ _ 

10^^^ eV m~^ s~^ sr~^ for both proton dominated com- 
positions and for ou r proton to iron transitio n scenario 
(see Figs. 4 and 8 of lElecerprit fc Allardll20Tlf ) . 

For single sources, iDecerprit fc AUardl (j2011| ) showed 
that the cosmogenic neutrino flux could be within 
reach of IceCube for powerful steady sources (see also 
lEssev et"ani2010f ). Only beamed sources (i.e., blazars) 
seem to satisfy the required luminosity condition (other- 
wise, the required power exceeds the Eddington power), 
but the neutri no flux is then dilut ed by the deflection of 
cosmic rays (iMurase "eFaIll20Tll ). In the case of tran- 
sient sources, the total received flux should be diluted 
by the ratio of the emission time to the spread in the 
arrival times due to the magnetic fields, Aig/St, which 
could lower the fiux of many orders of magnitude, pre- 
venting a ny detection. In the same toke n, as was dis- 
cussed in iGabici &: AharonianI (|2005D and IKotera et al.l 
()2011|) , the secondary gamma-ray emission (produced in 
the intergalactic medium) from a single transient source 
should be affected by dilution in time, and be below reach 
of next generation gamma-ray instruments. 

IMurase et all ()2009[) calculated that a promising 
amount of neutrinos would also be produced via hadronic 
interactions during the escape of ultrahigh energy pro- 
tons from the surrounding supernova envelope. These 
authors show the importance of muon and pion interac- 
tions on the baryonic envelope for the out-coming neu- 
trino flux, especially for interactions at the earliest times. 

Our current simulations do not take into account such 
interactions, hence an accurate evaluation of the neutrino 
flux associated to our scenario cannot be computed with 
the present tool. It can be noted however, that the over- 
all background neutrino flux that we would obtain would 
be sim ilar to the neutrino flux predicted bv IMurase et al.) 
()2009D for proton injection and at least about one order of 
magnitude smaller for iron injection. Our lower required 
source rate should not affect the neutrino flux, as its level 
is determined mainly by the energy injected in UHECRs 
above the pion production threshold energy. The su- 
pernova envelope opacity necessary to allow the escape 
of iron nuclei at the highest energies results indeed in a 
neutrino flux about o ne order of magnitude low er than in 
the case of protons (jMurase fc Beaconil[20T0l ). We also 
calculated that the primary injection of iron would lead 
to an enhanced flux (of a factor of a few) around PcV en- 
ergies because of the steeper overall UHECR spectrum 
generated after the escape. The level of neutrino flux 
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co uld still turn out to be fairly high (the flux predicted 
bv lMurase et al.|[2009l is significantly above the IceCube 
sensitivity), and leaves room for detection with the Ice- 
Cube experiment, either of single close-by sources born 
within ~ 5 Mpc, or of the diffuse background. A full cal- 
culation of the expected flux is needed to formulate more 
quantitative statements. 

The gamma-rays produced in the supernova envelope 
via hadronic interactions could cascade in turn on the 
same background and escape as photons in the TeV 
range. This process could possibly produce a bright tran- 
sient gamma-ray source, though the exact spectrum and 
its detectability have to be quantitatively calculated. A 
fraction of ultrahigh energy photons could also escape, 
that could be also observed as a transient source by ex- 
periments such as Auger or JEM-EUSO, for sources at a 
few megaparsec distances (|Muras i [20091) . A gain, these 
assertions need more careful investigations. 

Highly magnetized magnetars with fields B > 10^^ G 
should be strong emitters of gravitational waves. If pro- 
tons are injected by pulsars, the hard produced spectrum 
requires a specific distribution of pulsar parameters (of 
their initial rotation velocity and/or their magnetic field 
strength) to softe n the overall UHECR spectrum and fit 
the observations (jKoterair2011[ ). In s uch a ca s e, and for 
strong pulsar internal deformations, iKoteral (j2011[ ) ar- 
gues that a characteristic diffuse gravitational wave sig- 
nal would be produced, that could be detected with fu- 
ture generation detectors such as DECIGO or BBO. The 
present study shows however that these strong magnetic 
fields would induce a fast spin-down time that could not 
allow the escape of UHECRs in presence of a dense super- 
nova shell. The problem could be bypassed if the enve- 
lope is particularly under-dense, if particles could escape 
through a breech created by a proto-pulsar jet (but in- 
teractions with the radiative background would no longer 
be negligible in that case), or for e nvelope shredding sce- 
narios as invoked by lAronsI (j2003f ) . For the milder fields 
favored in our scenario {B ~ 10^^ G), the gravitational 
wave signal is expected to be lower by many orders of 
magnitude, far below the reach of any planned instru- 
ments. 

4.6. A Galactic scenario for UHECRs? 

We discuss in this section the scenario in which the 
major contributor to the currently observed UHECR flux 
are Galactic pulsars injecting iron, and possessing the pa- 
rameters required for iron escape at the highest energies. 
The pure proton injection case seems indeed difficult to 
reconcile with a continuous detection of UHECR events, 
given the short spread in their arrival times ^^cai at the 
highest energies and the lack of anisotropics toward the 
Galactic plane. The iron injection case could be more 
promising, in so far as Jtcai can be much longer than the 
detector exposure time, Texp, even at the highest ener- 
gies, for reasonable Galactic magnetic field strengths. 

We note Vg the birth rate of neutron stars in our Galaxy 
that satisfy the conditions for successful iron accelera- 
tion and escape at the highest energies (Eq. [57]). We 
recall that the 'normal' pulsar birth rate is of order 
J^Gai ~ 1/60 yr-i in the Galaxy, (jLorimeii 120081) . If the 



time interval between two births is shorter than the dis- 
persion of the arrival times i'^^ < StQ^i, then the flux 
of UHECRs should not depend on StQ^i and could be 
accounted for by the fraction ~ 10^^ — 10~^ of the pop- 
ulation of pulsars within our Galaxy. 

Now, if > Stoai, one may have zero (in which case 
the Galactic scenario does not stand) or only one source 
contributing to the observed Galactic UHECR flux. The 
UHECR fiux due to this source can be written: 

E'ME)=E^^j^^^ (29) 
~4x 103°eVV-2s-isr-i- 

^26^30.5*10 kpc 

This estimate implies that in this scenario, the cosmic- 
ray injection efficiency should be ^ 4 x 10~^ times lower 
not to overshoot the spectrum. 

In the scenario where extragalactic pulsars dominate 
the observed UHECR flux, ^z^s = /sJ^Gai, with ^ the 
iron injection efficiency. With /g < 0.01% as calcu- 
lated before, we would likely fall in the latter case, with 
> Stcai. In this scenario, there should be no Galac- 
tic source contributing currently, as otherwise, it would 
overshoot the observed spectrum. 

Note also that these flux estimates are subject to strong 
variations according to the structure and strength of the 
Galactic magnetic field. 

In both cases (single or many Galactic sources con- 
tributing), the energy spectrum should present a cut-off 
(at Edit) mimicking the GZK cut-off, due to the propa- 
gation of particles in the supernova envelope. The ankle 
feature would stem from the change in slope around the 
secondary proton peak. The overall spectral index could 
fit the observed one by a combination of the escaped 
spectrum and the propagation effects in the Galaxy. 

The chemical composition of UHECRs detected on 
Earth would slightly differ from the composition of par- 
ticles escaped from the supernova envelope. Protons 
would indeed disappear more quickly from the Galaxy 
than heavy elements around 10^^ eV, as Stcai of order of 
Texp at this energy. At lower energies, particles should 
still be able to remain confined, and a transition from 
light to heavy nuclei should still occur. 

Finally, the main weakness of this Galactic scenario 
lies in the expected anisotropy signature. A single source 
should lead to a noticeable spot of events in the sky at 
the highest energies, even for iron nuclei, unless the tur- 
bulent Galactic magnetic field is extremely strong. If 
many sources were contributing, they are also expected 
to trace the Galactic disk, but no such anisotropy has 
been observed in the UHECR data. 

5. CONCLUSIONS 

We studied the injection and escape of UHECRs from 
newly-born pulsars based on a Monte Carlo simulation 
of hadronic interactions and on a detailed examination 
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of the physical properties of supernovas envelopes. Our 
results show that protons and light elements at the high- 
est energies can traverse the envelope only for very di- 
lute envelopes. For pulsars embedded in supernovse with 
characteristics satisfying Eq. ([27]) . iron nuclei arc able to 
escape from the supernova envelope with energy above 
10^° eV. The escaped spectrum displays a transition 
from light to heavy composition at a few EeV, match- 
ing the recent Auger data. Due to the production of 
secondary nucleons, the escaped spectrum also presents 
a softer slope than the initial injected one, enabling a 
good fit to the observations. 

Under the assumption that unipolar induction accel- 
eration can take place in the neutron star winds, two 
conditions ensure the compelling adequacy of the sce- 
nario of production of UHECR by neutron stars with 
the observed data: 



1. that a fraction fs < 0.01% of extragalactic super- 
nova: give birth to pulsars with sub-millisecond 
periods, and dipole magnetic field in the range 

1012-13 G, 

2. that a successful injection of heavy nuclei is occur- 
ring at the acceleration site of these objects. 

We discussed that this double condition can be rea- 
sonably fulfilled. Indeed, about 2% of young 'normal' 
pulsars arc inferred to have initial parameters close to 
fulfilling condition 1. The low value of fs also permits 
poorer injection efficiencies and compensates for the nar- 
rowness of the allowed parameter range. Condition 2 is 
naturally favored in neutron stars that have heavy nuclei 
rich surfaces. Should these two conditions be fulfilled, 
the main UHECR observables, namely, the energy spec- 
tra, composition, and arrival directions, would be consis- 



tent with the lates t Auger data ()Abraham et al.ir2010bl lal: 
lAbreu et al.l[2010l) . 

If criterion 2 is not met, and only protons or light el- 
ements can be accelerated, the neutron star scenario is 
viable only for dilute surrounding envelopes, or if mech- 
anisms shredding or piercing the envelope are at play to 
enable the escape of particles. 

In the iron injection case, the birth of such an object 
within our Galaxy would be noticeable in the number 
of detected events, only for very close-by sources (at 
^ 2kpc). Such a source could lead to a distinct increase 
of the observed UHECR events starting (Jtcai ^ 70yrs 
x(Z/2kpc) after the birth (for a source located at with 
parameters chosen in Eq. [^5|) . and that would last for 
Jicai- If a pulsar birth were observed today, the proton 
injection case and/or the production of neutrons by in- 
teractions in the direct environment of the source would 
lead to a significant burst of UHECR events in the sky, a 
fraction of year later. The birth rate of neutron stars sat- 
isfying our criteria inside our Galaxy is however expected 
to be as low as ~ 5 x 10"'' yr-^. 

Other signatures can be expected, such as a non- 
recovery of the energy spectrum above £^cut,Fo ~ 
2Q20.5 Qj. ^Yjc precise measurement of the cosmic ray 
composition at high energies. Large exposure instru- 
ments such as Auger North or JEM-EUSO would allow 
to make such measurements and probe this scenario. 
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APPENDIX 

ALTERNATIVE MECHANISMS TO ESCAPE THE SUPERNOVA ENVELOPE? 

Recent works have shown the possibility that the confining pressure of the toroidal magnetic field could collimate 
the proto-magnetar wind along its polar axis, and drive a jet th at has the properties of long gamma-ray bursts jets 

omissarov fc Barkovll2007l : iBucciantini eTall [20071 [200l l2009l ). This scenario opens up the possibility that cosmic 
rays be accelerated via magnetic reconnection or Fermi acceleration inside the proto-magnetar jet, and escape through 
the pierced superno va envelope. The cas e of nuclei escaping through a je t has been discussed s emi-analytically in the 
context of GRBs bv lMurase et al.l (|2008f ) and for proto-magnetar jets bv lMetzger et al.l ()2011aD . 

However, mildly mag netized pulsars could not have the collimation power to produce a jet. As discussed in 
IBucciantini et ahl (|2007[ ). the collimation becomes significant for values of the ratio of the Poynting flux to the to- 
tal energy at the termination shock of the wind, £'mag/£'tot ^ 0.2, at times i ^ 10 — 100 s. The conversion of magnetic 
energy into kinetic energy in relativistic outflows at large radii is uncertain. For the mild magnetic fields and high ro- 
tation velocities that we consider, the magnetization at the light cylinder reads (Tl = 4:iJ?il^ / (Mc^) ~ 20/i|Q ^^lf/Mg,5, 
where the mass loss rate at i ^ 30 s, Afg.5 = M /10~^'^ Mq, was inferred from Appendix A2 of iMetzger et al.l ()2011bD . 
The value of ctl should increase steeply with the fast mass loss rate around t ~ 60 s and later (when the wind becomes 
transparent to neutrinos). For such high ctl 3> 1, it is plausible that magnetic dissipation in the relativistic outfiow 
out to the termination s hock leads to a low £'mag/£'tot at these distances (jCoronitil 119901 : iLvubarskv &: Kir3 120011 : 
iKirk fc Skiaeraas"eiill2003f l. not allowing the formation of a jet. Studies of the Crab Pulsar wind nebula show indeed 
that i^mag/^^tot 10"2 at large radii (jKennel fc CoronitillltsllBegelman fc Lilll992D . 

One may note that, even in a scenario where a jet were produced (for more strongly magnetized neutron stars). 
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our conclusions would still apply, for particles that would not be injected in the direction of the jet, but in the other 
sectors. Such particles would have to cross the expanding supernova shell, and would have the same fate as in the 
present framework. 

Ano ther mechanism to bypass the problem of UHECRs crossing the dense supernova envelope was invoked bv I AronsI 
(|2003[) . who proposed that the supernova envelope be disrupted by the magnetar wind. Such phenomena have never 
been observed, neither in magnetar envelopes, nor in rotation-powered pulsar envelopes. 

Finally it is also pos sible that millisecond pulsars can be born in Accretion- Induced Collapse of white dwarfs 
(|Frver et all 1 1991 [20091) . The escape scenario is different and will be studied in our future work. 
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